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摘 要 :为 了 探究 单 层 柱 面 网 壳 结 构 受 冲击 作用 下 的 动态 响应 ,验证 有 限 质点 法 模拟 网 过 结构 冲击 
动态 响应 的 准确 性 ,建立 了 7 个 冲击 工 况 的 有 限 质 点 法 模型 ， 采用 以 对 称 罚 函 数 为 基础 的 双 球 碰撞 
模型 为 冲击 模型 ,根据 能 量 守恒 定律 确定 冲击 模型 中 碰撞 刚度 系数 ,引入 Cowper-Symonds 本 构 模 
型 ,获取 结构 动力 响应 ,并 将 模拟 结果 与 试验 结果 对 比 。 结 果 表 明 :1) 有 限 质点 法 得 到 的 冲击 力 时 
程 变化 曲线 与 试验 结果 相近 ;2 ) 由 于 未 考虑 结构 损伤 因素 ,有 限 质 点 法 所 得 节点 位 移 略 小 于 试验 
值 ,但 所 获得 的 结构 失 稳 动量 与 试验 一 致 ;3) 冲击 产生 的 网 过 应 变 能 更 易 沿 着 网 党 斜 杆 传递 ,同时 
冲击 区 杆 件 ,刚度 较 大 区 域 如 支 座 附近 的 杆 件 为 结构 的 主要 受 力 构件 。 有 限 质 点 法 可 以 较为 准确 
的 判断 单 层 柱 面 网 这 结构 受 冲击 作用 的 动态 响应 ,可 为 提高 相关 结构 抗 冲击 性 能 提供 理论 依据 。 
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Abstract: To study the dynamic response of single-layer cylindrical latticed shell under impact and 
examine the accuracy of finite particle method( FPM) in simulating the dynamic response of latticed shell 
structure ,7 FPM models were established based on 7 impact conditions , which were adopted contact model 


of double ball deriving from symmetrical penalty function and were introduced Cowper-Symonds 
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constitutive model. The law of energy conservation was used for obtaining impact stiffness coefficients in 


contact model. The results of simulation were used to reveal the dynamic response of structure and were 


compared with experimental ones. The results show that 1 ) time history curves of impact force were close 


those of experimental ones,2) node displacements got in FPM were smaller than those of experiment due 


to failure of considering damage factors of latticed shell ,but its destabilizing momentum was the same as 


that of experiment ,3 ) strain energy which was generated from impact had a tendency to transfer through the 


braces. Meanwhile ,members in impact area and high stifftness area, such as adjacent supports area, were 


the main forced components. In total ,FPM is an accurate method to simulate dynamic response of single- 


layer cylindrical latticed shell under impact, which can provide theoretical basis for improving impact 


resistance of similar structures. 


Key words:impact resistance of single-layer cylindrical latticed shell ;finite particle method ;test of single- 


layer cylindrical latticed shell; dynamic response;comparative analysis 


空间 网 格 结构 具有 跨度 大 、 造 型 优美 .质量 轻 、 
便于 标准 化 等 特点 ,在 机 场 、 航 站 楼 、 体 育 场馆 等 公 
共 建 筑 中 大 量 使 用 。 为 了 保证 结构 服役 阶段 的 安全 
性 ,许多 学 者 对 网 格 结构 在 重力 荷载 `. 风 和 荷 载 `. 地 震 
作用 下 的 力学 性 能 及 抗 连续 性 倒塌 性 能 进行 研究 。 
相 较 于 承载 力 不 足 .地震 作用 等 因素 造成 的 空间 网 
格 结构 损伤 事故 ,冲击 损伤 事故 发 生 的 概率 较 低 , 但 
后 果 往 往 较 为 严重 。 如 1980 年 ,一 架 4 座 飞机 险 入 
巴西 圣保罗 展览 中 心 的 网 架 屋顶 ,造成 严重 的 财产 
损失 号 。 研 究 空间 网 格 结构 受 冲 击 的 动态 力学 特 
性 ,分 析 其 变形 机 理 具 有 十 分 重要 的 意义 。 

空间 网 格 结构 抗 冲击 性 能 的 数值 模拟 研究 主要 
基于 有 限 元 法 , 常 采用 拉 格 朗 日 乘 子 法 和 对 称 罚 画 
数 法 冲击 模型 ,前 者 可 用 于 隐 式 求解 ,后 者 在 显 式 和 
隐 式 求解 中 均 可 使 用 2 ,二 者 都 可 在 有 限 元 软件 平 
台 实 现 。 文 献 [3] 和 文献 [4] 均 用 ANSYS/LS-DYNA 
研究 了 网 格 结构 受 冲击 时 的 动态 响应 。 但 为 保证 收 
敛 ,有 限 元 法 对 时 间 步 和 增 量 步 的 要 求 较为 严格 中 。 

2004 年 TING 等 1, 提出 了 基于 向 量 式 结构 与 
固体 力学 的 数值 模拟 方法 。 喻 莹 ”在 此 基础 上 面向 
结构 工程 领域 提出 了 有 限 质点 法 (finite particle 
method,FPM) 。 相 较 于 有 限 元 法 ,有 限 质点 法 在 分 
析 接 触 碰撞 问题 时 不 需要 迭代 求解 ,在 每 一 时 间 步 
内 单独 计算 接触 力 和 接触 内 力 ,并 反作用 于 质点 ,无 
需 修 正 刚度 方程 ,并 以 显 式 中 央 差 分 法 作为 求解 策 
略 ” 。 喻 莹 等 用 有 限 质 点 法 ,采用 对 称 罚 函 数 冲 
击 模 型 模拟 了 柔性 体 和 刚性 体 边界 碰撞 时 的 形态 变 
化 ,以 及 空间 梁 磁 撞 时 的 形态 变化 ,计算 结果 较为 准 
确 。 有 限 质点 法 亦 用 于 网 格 结构 受 冲 击 性 能 的 数值 
分 析 , 喻 莹 等 \ 采用 理想 弹 塑 性 本 构 模 型 和 Mises 


屈服 准则 ,用 有 限 质 点 法 模拟 K8 网 壳 受 冲击 后 的 连 
续 性 倒塌 过 程 。 

但 目前 基于 有 限 质点 法 研究 网 格 结构 受 冲 击 作 
用 的 动态 响应 与 试验 数据 对 比较 少 ,同时 所 采用 的 
材料 本 构 模 型 基于 材料 的 静态 力学 性 能 ,而 冲击 问 
题 需 考虑 材料 在 冲击 荷载 作用 下 届 服 强度 和 瞬时 应 
力 随 应 变 率 的 变化 "1。 为 此 ,本 研究 将 考虑 动态 屈 
服 强度 和 应 变 率 关系 的 本 构 模 型 引入 单 层 柱 面 网 壳 
的 有 限 质 点 法 数值 分 析 , 主要 研究 内 容 为 :1 ) 将 基于 
对 称 罚 函数 法 的 冲击 模型 和 适用 于 金属 冲击 的 
Cowper-Symonds 本 构 模 型 编 人 有 限 质点 法 程序 ;2 ) 
建立 单 层 柱 面 网 过 结构 有 限 质 点 法 模型 ,研究 网 壳 
结构 受 冲击 作用 的 动态 力学 响应 ,并 与 试验 结果 对 
比 ;3) 从 应 用 角度 ,研究 单 层 柱 面 网 过 在 冲击 荷载 作 
用 下 的 能 量 传递 路 径 ,探究 结构 传 力 特性 , 找 出 关键 
受 力 构件 ,为 类 似 结 构 的 抗 冲击 研究 和 工程 应 用 提 
供 依据 。 


1 有 限 质点 法 的 基本 理论 


有 限 质 点 法 将 结构 离散 成 一 组 有 限 数目 的 质 
点 ,以 质点 为 最 小 单元 描述 结构 的 性 质 ,单元 仅 起 到 
连接 节点 的 作用 。 该 方法 以 牛顿 第 二 定律 为 控制 方 
程 ,建立 关于 时 间 的 常 微分 方程 式 ,将 静 力 和 动力 问 
题 中 质点 的 受 力 状态 视 为 准 动 平衡 状态 和 动 平衡 状 
态 , 其 中 毅力 问题 通过 设置 虚拟 阻尼 系数 实现 。 式 
(1) 和 式 (2) 为 具有 6 个 自由 度 的 空间 质点 / 的 控 
制 方程 。 


mx =F,= 下 和 +F (1) 
1,0,=M,= M” +M” (2) 
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式 中 :mm) 为 点 了 的 质量 ,包括 点 /自重 和 相连 单元 所 
提供 的 等 效 质量 ;*, 和 0, 为 点 / 的 线 加 速度 向 量 和 
角 加 速度 向 量 ; 卫 为 点 的 质量 惯性 向 量 ;为 点 J 
的 合力 向 量 ,包括 作用 在 点 /的 集中 力 向 量 了 PY 与 
点 了 相连 单元 所 提供 的 等 效 内 力 向 量 Fy 等 ;MM 为 
点 的 合 弯 矩 向 量 , 包 括 作 用 在 点 7 的 外 力矩 向 量 
M7 与 点 相连 单元 所 提供 的 等 效 弯 矩 向 量 
MY” 等 

变形 机 制 是 有 限 质 点 法 的 核心 ,有限 质 点 法 
采用 虚拟 逆向 运动 和 正 向 运动 ”的 方法 得 到 固体 
结构 扣除 刚体 位 移 后 的 纯 变 形 , 并 由 此 计算 作用 在 
质点 的 上 内 力 和 在 单元 中 产生 的 应 力 和 应 变 , 而 内 
力 和 位 移 之 间 的 关系 则 通过 单元 材料 本 构 方 程 
建立 。 

得 到 内 力 后 ,有 限 质 点 法 依据 显 式 中 央 差 分 法 
计算 下 一 时 刻 质点 的 位 移 , 见 式 (3) 和 式 (4) 。 如 此 
往复 ,直至 时 间 终 止 。 


下 ， 和 
Xnr1 = CI ne +2C — CX, 1 (3) 


M, 及 
0 = (Ci 7 +2C10, - C0,_1 (4) 


式 中 :为 时刻 作 用 在 质点 上 的 合力 向 量 和 合 弯 
和 矩 向 量 ;下 标 n -1,n 和 n+1 分 别 代 表 i, -1、t 和 
t, +1 时刻;h 为 时 间 步 长 ;Cl 和 C, 分 别 为 虚拟 阻尼 
系数 。 当 =0 时 ,zx ,和 60. 没有 实际 的 物理 含义 ， 
因此 根据 式 (5) 和 式 (6) 来 计算 x, 和 0。 
7 . Fo +tf, 


Xi =X0 十 帮工 十 三 (5) 
? M+ 

EY (6) 
2 到 


式 中 ,f 和 m, 为 设置 的 虚拟 阻尼 力 和 虚拟 质量 。 
2 冲击 试验 


本 课题 组 刘 建 伟 ' 进行 了 冲击 作用 下 单 层 柱 面 
网 壳 动 态 响 应 的 试验 ,以 研究 其 动态 响应 和 动力 稳 
定性 。 

试验 选用 规格 为 B10 x2 的 无 缝 钢管 ,弹性 模 量 
=1.87 x10”MPa, 剪 切 模 量 C =7.19 x10” MPa, 局 
服 强 度 为 262 MPa, 极 限 强 度 为 368 MPa, 泊 松 比 为 
0.3, 单 位 长 度 质量 mw =0. 222 kg/m。 杆 件 通 过 尺寸 
为 100 mm x 100 mm x8 mm 的 节点 板 焊接 在 一 起 。 
网 壳 以 贸 接 的 方式 与 反 力 架 相 连 。 冲 击 物 为 质量 分 
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别 是 20、5 2.18 kg 的 硬 质 圆 钢 。 网 党 跨 度 为 1 500 mm， 
矢 高 为 375 mm ,纵向 长 度 为 2 000 mm ,在 点 妃 放 置 
YZC-526/1 000 kg 应 变 式 力 传感器 ,通过 应 变 仪 采 
集 冲 击 力 , 并 在 点 4、B、C、D、E、F 放置 DH808 拉线 
式 位移 传 感 咒 。 图 1 和 图 2 分 别 为 现场 装置 图 和 试 
验 装置 示意 网 。 

力 传感器 


试验 样品 


位 移 传感器 
图 1 试验 现场 装置 图 
Fig.1 Diagram of field test device 
冲击 物 通过 设置 于 反 力 架 上 横梁 的 滑轮 装置 ， 
以 自由 落体 的 方式 冲击 网 壳 节 点 B, 如 图 2(b) 所 
示 。 通 过 调整 冲击 物 的 下 落 高 度 来 改变 冲击 速度 ， 
共 包 含 表 1 所 示 的 7 种 工 况 。 其 中 冲击 速度 是 冲击 
物 和 网 过 接触 时 的 瞬时 速度 ;WC 代表 工 况 。 


1 500 


KS 


2 000 | 


控制 台 


(b) 冲击 力 加 载 装 置 示意 图 


图 2 试验 加 载 示 意图 
Fig.2 Test loading scheme 
考虑 到 加 工时 可 能 出 现 的 尺寸 误差 ,试验 前 对 
结构 杆 件 实际 长 度 进行 3 次 测量 校正 ,测量 结果 取 
平均 值 , 作 为 数值 模拟 的 模型 几何 参数 。 
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表 1 各 工 况 计算 参数 


Tab. 1 Experimental conditions 


冲击 点 。 WC 冲击 物质 量 / 冲击 物 速度 / 冲击 动量 
kg (ms ) (kg ms  ) 
B 1 2.18 2 4.36 
B 2 2.18 2.844 6.2 
B 3 2.18 4 8.72 
B 4 5 1 5 
B 5 a 2 10 
B 6 5 4 20 
B 7 20 6 120 


3 有 限 质点 法 程序 编制 


为 简化 分 析 ,假设 冲击 物 为 刚性 体 ;冲击 过 程 产 
生 的 热能 和 摩擦 忽略 不 计 ; 发 生 多 次 冲击 时 ,第 一 次 
冲击 的 结构 动态 响应 最 为 关键 , 故 本 人 研究 只 考虑 一 
次 冲击 。 


3.1 建立 冲击 模型 


基于 本 课题 组 此 前 研究 ,采用 以 对 称 罚 函数 
为 基础 的 双 球 碰撞 模型 。 

对 称 罚 函数 原理 为 ,每 一 个 时 间 步 都 要 检查 冲 
击 物 和 被 冲击 物 之 间 是 否 发 生 穿 透 ,没有 则 不 作 任 
何 处 理 , 认 为 冲击 力 为 零 ; 如 果 发 生 穿 透 , 则 假想 冲 
击 物 和 被 冲击 物 之 间 有 一 个 弹簧 (图 3) 阻 碍 穿 透 发 
生 。 假设 弹 得 的 刚度 为 ,可 建立 冲击 力 计算 公 
式 , 即 


人 ，A<0 (7) 


| 3 
Ce 


图 3 双 球 碰撞 模型 
Fig.3 Contact model of double ball 
图 中 :5 为 假想 弹 复 的 刚度 系数 , 即 碰撞 刚度 系数 ; 
0 为 被 冲击 物 ;0, 为 冲击 物 ;di、ds 分 别 为 被 冲击 物 
和 冲击 物 的 质点 到 接触 面 的 距离 ;A 是 两 球 外 表面 
间 的 距离 ,车 A=0 表明 没有 发 生 冲 击 行为 ,反之 则 
发 生 冲 击 行为 。 
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3.2 确定 碰撞 刚度 系数 


根据 式 (7) ,确定 碰撞 刚度 系数 是 计算 冲击 力 的 前 
提 。 本 节 从 能 量 角度 给 出 确定 碰撞 刚度 系数 的 方法 。 
冲击 过 程 的 能 量变 化 可 分 为 2 个 阶段 "i。 第 
1 阶段 时 间 为 毫秒 级 ,冲击 物 的 动能 部 分 转化 成 网 
过 冲击 区 的 动能 ,此 时 非 冲 击 区 的 动能 和 和 结构 应 
变 能 较 小 ;第 2 阶段 中 冲击 区 的 动能 迅速 向 非 冲 击 
区 传递 ,结构 应 变 能 逐渐 增加 ,最 后 结构 的 动能 通过 
结构 的 势能 和 应 变 能 耗 散 ,这 一 过 程 的 时 间 以 秒 计 。 
针对 高 速 冲击 ,两 阶段 能 量变 化 特征 较为 明显 ,可 以 
第 1 阶段 未 为 临界 点 将 网 壳 整 体 动 量 和 能 量 守 恒 转 
化 为 冲击 区 的 动量 和 能 量 守 恒 。 但 对 于 低速 冲击 ， 
非 冲 击 区 在 第 1 阶段 就 已 产生 明显 的 动力 响应 5 ， 
上 述 方 法 则 不 再 适用 , 故 采 取 式 (8) ,根据 整体 结构 
的 能 量 守恒 求解 碰撞 刚度 系数 , 即 
A 取 +A 有 =A 了 +AE +AE, (8) 
式 中 :Acev 和 Akcw 分 别 为 冲击 物 势能 改变 量 和 动 
能 改变 量 ; AW.s、AERs 和 AE; 为 分 别 为 结构 势能 改 
变量 动能 改变 量 和 应 变 能 改变 量 ; 式 (8) 适用 于 未 
发 生 构件 断裂 的 结构 。 对 每 一 个 冲击 工 况 , 均 有 唯 
一 的 结构 碰撞 刚度 系数 与 之 对 应 。 硅 数值 分 析 
所 输入 的 碰撞 刚度 系数 .不 等 于 4, 则 式 (8) 不 再 
守恒 , 令 等 式 左边 与 右边 的 差 值 为 冲击 前 后 系统 能 
量 差 AWW。 通 过 数据 拟 合 ,以 A 勾 =0 时 所 对 应 的 天 
为 结构 碰撞 刚度 系数 〖,。 确 定 磁 撞 刚 度 系 数 的 流 
程 如 图 4 所 示 。 


输入 Kk 


AswAHewAFasAFsAPa 


| 


计算 A 矿 ， 并 输出 A 矿 和 大 
获取 多 组 数据 
拟 合 AW-k. 


图 4 计算 碰撞 刚度 系数 流程 图 


Fig.4 Flow chart of calculating impact stiffness efficient 
3.3 ”本 构 模 型 的 选取 


相 较 于 静 力 学 问题 ,冲击 问题 需 考虑 惯性 效应 
和 应 变 率 效 应 " ,常用 的 金属 动态 力学 行为 本 构 模 
型 包括 Cowper-Symonds 模型 、Johnson-Cook 模型 、 
Bammann 模型 和 Zerilli-Armstrong 模型 。Cowper- 
Symonds 模型 建立 了 应 变 率 和 材料 动态 屈服 强度 的 
关系 。Johnson-Cook 模型 在 此 基础 之 上 考虑 了 温度 
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对 动态 届 服 强度 的 影响 。Bammann 模型 的 基本 方 
程 需 要 包括 损伤 参数 .内 变量 等 在 内 的 多 种 材料 参 
数 吕 ] 。Zerilli-Armstrong 模型 基于 位 错 动 力学 ,考虑 
了 温度 的 影响 ,适用 于 体 心 立方 模型 和 面 心 立方 模 
型 中 中。 鉴于 数值 模拟 参照 试验 在 室温 环境 下 进行 ， 
故 不 考虑 温度 对 冲击 试验 的 有 影响。 同时 Bammann 
模型 所 需 材料 参数 较 多 , 故 选 择 Cowper-Symonds 模 
型 ,该 模型 在 网 壳 结 构 受 冲击 研究 中 被 广泛 应 用 , 模 
型 公式 为 


-= 


ou = (oo+E,s,)|1+ 加 ] (9) 


式 中 :ou 和 oo 分 别 为 动静 态 屈服 应 力 ;2 和 8 分 别 为 
等 效 塑 性 应 变 和 等 效应 变 率 ;为 塑性 硬化 模 量 ;C 
和 p 为 常数 。 针 对 参照 试验 ,C 和 pp 分别 取 40 
和 553]。 

冲击 力学 模型 程序 编写 思路 参照 文献 [ 19 ] , 根 
据 式 (10) 判 断 材料 是 否 进 入 塑性 状态 。 若 B >0， 
则 进入 塑性 状态 ,反之 则 为 弹性 状态 。 在 判断 前 , 认 
为 单元 统一 处 于 弹性 阶段 ,应 力 应 变 根据 弹性 理论 
计算 ,此 时 所 得 应 力 为 预测 应 力 。 当 满足 B >0 时 ， 
单元 进入 塑性 阶段 ,此 时 需 根 据 式 (11) 更 新 单元 应 
力 应 变 '” ,使 之 满足 塑性 要 求 ,更 新 后 的 应 力 称 为 
真实 应 力 。 


DD=o0-0o (10) 
站: 
AT (11) 


OrAt 

式 中 :下 标 t+At 代表 时 刻 ;为 等 效应 力 ;为 动态 

屈服 强度 ;yw 和 cpuosw 为 真实 应 力 和 预测 应 力 ; 

C 为 剪 切 模 量 ;S;j ,ws 为 偏 应 力 分 量 ;AA 为 塑性 流 

动因 子 增 量 。 各 参数 的 详细 计算 方法 见 文 献 [19]。 
结构 材料 采用 分 段 线性 弹 塑性 模型 , 见 图 5。 


40 


Oij,t+At = Oj,tr,t + At -3GAA 


0 0.05 0.10 0.15 0.20 025 
应 变 


图 5 材料 分 段 线性 弹 塑性 模型 


Fig.5 Piecewise linear elastoplastic model of material 


第 40 卷 


3.4 程序 编制 


程序 编制 流程 如 图 6 所 示 。i% 时 刻 模 型 包含 网 
这 质点 坐标 和 冲击 物 坐 标 。 结 构 在 加 时 刻 的 三 维 模 
型 如 图 7 所 示 。 


输入 太 、# 时 刻 静 载 和 单元 内 力 


根据 式 (7) 计 算 
a 


击 何 蕴 


图 6 程序 编制 流程 图 
Fig.6 Flow chart of programming 


图 7 有 限 质 点 法 三 维 模型 
Fig.7 3D model of FPM 


4 结果 与 分 析 


4.1 确定 碰撞 刚度 系数 


本 节 基 于 有 限 质点 法 利用 MATLAB R2019a 数 
学 分 析 软 件 模拟 上 述 试验 过 程 ,并 对 比 模拟 结果 。 
文献 [9] 建 议 取 1 x10 ~1 x10™s 作为 碰撞 问题 
的 时 间 步 长 ,本 研究 取 1 x107s, 模 拟 试验 中 7 个 
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工 况 。 
以 工 况 5 和 工 况 6 为 例 展现 碰撞 刚度 系数 的 求 
解 过 程 , 表 2 记录 了 第 一 阶段 末 工 况 5 和 工 况 6 选 
取 不 同 磁 撞 刚 度 系数 所 获取 的 能 量 参 数 , 建 立 AW-k. 
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拟 合 曲 线 , 确 定 结 构 碰 撞 刚度 系数 4。 根据 图 8 拟 
合 结果 , 工 况 5 和 工 况 6 的 结构 碰撞 刚度 系数 分 别 
为 90 934.241 86 121 318. 128 9 Nm ,其余 工 况 的 求 
解 结果 见 表 3。 


表 2 工 况 5 和 工 况 6 不 同 碰撞 刚度 系数 所 对 应 的 能 量 参数 


Tab.2 System energy corresponding to different impact stiffness coefficients in case $5 and case 6 


WC k./( N.m-! ) AWcw + AEcw/J AEi/] AEs/J AWcs/J AW/J 
5 80 000 -3.3992 0.722 4 2.5640 0.0606 -0.0522 
90 000 -3.8109 1.1100 2.6242 0.0763 —0.0004 
100 000 -4.1632 1.330 5 2.7818 0.0807 0.029 8 
110 000 —4.4549 1.6270 2.8304 0.0850 0.0875 
120 000 -4.6927 2.0000 2.768 5 0.0889 0.1647 
6 100 000 -14.913 0 2.6027 11.9100 0.0350 -0.3660 
110 000 -16.2190 2.8706 13.0750 0.0458 -0.2280 
120 000 -17.3640 3.288 1 13.9800 0.0560 -0.0400 
130 000 -18.3870 4.1117 14.363 0 0.0661 0.1540 
140 000 -19.1740 3.999 5 15.279 0 0.058 4 0.1630 
02 25/20.575 =0. 122 <0.25, 故 有 限 质 点 法 得 到 的 半 
01 正弦 波 对 结构 动态 响应 的 影响 不 大 。 
0.2 
0 
0 Ea 
S§ -0.1 -0.2 
0.2 到 0.4 放 
= Se a a 
和 WC5 模 拟 结果 A 一 一 工 况 1(FPM) 
WC5 拟 合 曲线 = -0.6 
-03 A Wc6 模 拟 结果 革 MN ) en) 
一 一 Wc6 拟 合 曲线 os. \\ 属 一 一 工 况 3FPM) 
RY 1 ee ly 
80 000 100 000 120 000 140 000 -站 |r - 一 - 工 况 2( 试 验 ) 
Kk/(N*-m') - 一 - 工 况 3( 试 验 ) 


图 8 工 况 5 和 工 况 6 的 AW-h. 的 拟 合 曲 线 
Fig.8 Fitting curve of collision stiffness coefficient and 


system energy in condition 5 and 6 
冲击 力 时 程 变 


如 表 3 所 示 , 从 冲击 力 幅 值 角度 看 ,有 限 质 点 法 
和 试验 误差 范围 在 -20.29% ~14.98% 之 间 , 且 有 限 
质点 法 得 到 的 冲击 力 幅 值 和 动量 关系 与 试验 一 致 。 

图 9 是 有 限 质点 法 和 试验 得 到 的 以 工 况 1 ~3 
冲击 力 时 程 曲 线 对 比 图 。 由 对 比 结果 可 知 , 有 限 质 
点 法 得 到 的 时 程 曲 线 近似 半 正 弦 曲 线 , 试 验 曲线 近 
似 三 角形 。 不 同 的 冲击 荷载 形式 会 对 结构 反应 产生 
不 同 的 影响 。 然 而 , 当 荷 载 脉冲 长 度 比 w7, 即 冲击 
力作 用 时 间 和 结构 固有 自 振 周 期 的 比值 小 于 0. 25 
时 ,结构 动态 响应 主要 依赖 于 作用 冲 量 , 冲 击 荷 载 形 
式 对 其 影响 不 大 。 根 据 试验 ,3 种 工 况 冲击 持续 最 
长 时 间 为 25 ms, 且 用 3D3S 14 计算 结构 的 固有 周期 
为 205. 75 ms, 由 此 可 得 结构 荷载 脉冲 长 度 比 为 


4.2 


0 10 20 30 40 50 60 70 
时 间 /ms 
图 9 工 况 1~3 冲击 力 时 程 曲线 对 比 
Fig.9 Comparison of time history curves of 


impact force under conditions 1 ~3 


表 3 不 同 工 况 的 冲击 力 幅 值 


Tab.3 Amplitude of impact force under different conditions 
kx/ 冲击 力 幅 值 
WC 误差 /% 
(N.m-!) FPM/N 试验 值 /CN 

1 25 064.346 28 448.3682 562.5 一 20.29 
纪 30 238.684 02 675.3461 700.7 -3.62 
3 34 151.384 87 983.0852 1 139.6 =13573 
4 9 302.052 86 262.7988 254.6 3.22 
El 90 934.241 86 1041.007 905.4 14.98 
6 121 318. 1289 2 098.4332 2 025 3.63 
7 3 996.462 3 497 14.28 


28 820. 820 82 


结构 动力 失 稳 分 析 
结构 动力 失 稳 分 析 分 为 2 个 方面 ,其 一 为 同一 
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工 况 下 结构 不 同位 置 的 位 移 响应 ,其 二 为 综合 多 种 。” 络 发 生 大 幅度 增加 ,说 明 结 构 在 动量 处 于 20 ~ 

工 况 ,得 到 结构 的 临界 失 稳 动量 。 120 kg' m/s 时 发 生动 力 失 稳 ;2) 不 同 冲 击 动量 下 的 
图 10 为 有 限 质点 法 和 试验 得 到 的 工 况 7 点 C、 位 移 包 络 以 5 kg: m/s 最 小 ,以 120 kg m/s 最 大 ,中 


DEF 的 位 移 时 程 曲线 。 结 果 表 明 , 有 限 质 点 法 得 间 从 小 到 大 排列 依次 为 4.36 .6.2、10 和 20 kg:m/s; 
到 的 点 C 和 点 D 的 竖 向 位 移 数 值 偏 小 ,但 节点 位 移 3) 数 值 模拟 好 试验 均 表明 冲击 动量 为 5 kg . m/s 时 
变化 呈现 出 相同 的 规律 :1) 点 C 的 竖 向 位 移 最 和 完 达 ”的 结构 动力 响应 要 小 于 4.36 kg . m/s 下 结构 的 动 
到 最 大 值 ,其 次 是 点 D、 点 和 点 ;2) 在 冲击 作用 力 响应 ,这 是 表明 虽然 前 者 的 冲击 物质 量 较 大 ,但 冲 
下 ,4 个 节点 的 位 移 时 程 曲线 呈现 出 明显 的 四 陷 , 位 击 速 度 较 小 ,在 一 定 程度 下 ,速度 对 结构 动力 响应 的 
移 迅速 增 大 ,冲击 结束 后 ,节点 在 平衡 位 置 处 上 下 振 更 为 显著 。 


动 ;3) 在 结构 纵向 ,点 C 的 位 移 响 应 要 大 于 点 D 的 20 
位 移 响 应 ,但 由 于 支 座 约束 的 影响 ,在 跨度 方向 , 距 y 
离 冲击 点 较 近 的 点 的 位 移 响应 要 小 于 与 冲击 点 距 , 
a 目 
离 较 远 的 点 了 的 位 移 啊 应 。 上 
40 深 -10 
芭 
区 了 一 一 大 4.36 kg * m/s 
如 0 和 ZL —— [EE5 kg°* m/s 
加 0 所 天 一 一 大 62kg * m/s 
后 一 一 大 10kg * m/s 
汉 _40 = i kg * m/s 
三 —— [E120 kg * m/s 
| -50 
了 0 2 50 75 I 125 T1530 73 200 
想 时 间 /ms 
(a) 点 C(FPM) 
10 
0 50 100 150 200 250 300 3 
时 间 /ms 0 
(a) FPM 
B=5 
30 一 一 点 C Re 
一 一 点 D 广 _10 
2_ 一 点 有 生 时 —— [4.36 kg * m/s 
10 一 一 点 湖 对 —— [5 kg * m/s 
8 or 站 -20 一 一 大 6.2 kg * m/s 
这 0 2 2 —— E10 kg°* m/s 
革 _j0 | 二 一 一 大 20 kg * m/s 
= ES vy, 证 ——El20kg*m/ 
=90 2 DA | | ff 30 g* m/s 
-30 /Ey \ fp 0 2 50 75 100 125 150 175 200 
10 Sa CN WO 时 间 /ms 
Y Vv (b) 点 D(FPM) 
0 50 100 150 200 250 300 30 
时 间 /ms 
(b) 试验 20 
10 “ 工 况 7 节点 位 移 时 程 曲 线 对 比 We 
Fig. 10 Comparison of time history curves of 上 
displacement under condition 7 认 
文献 [20] 给 出 判断 单 层 网 壳 结 构 动 力 失 稳 的 2 人 
个 指标 :第 1 指标 通过 判断 以 网 壳 结构 运动 的 位 移 用 sem 
、 、 一 一 三 6.2 kg* m/s 
包 络 面 上 的 最 大 位 移 来 确定 失 稳 , 当 动 量 达 到 某 一 -30 * —— [10 kg . m/s 
4 i i i A 有 一 一 大 20 kg * m/s 
数值 时 ,结构 位 移 包 络 线 会 突然 增 大 ;第 2 指标 通过 -40 一 一 户 120 kg * m/s 
判断 特征 节点 竖 向 位 移 来 确定 失 稳 。 通 过 特征 节点 0 2 ISO 75200 
时 间 / 
竖 向 位 移 随 动量 的 变化 反映 网 壳 结 构 由 局 部 失 稳 到 四 0 
iv 时 2 =| vy 2 三 > 妈 
整体 失 稳 的 过 程 。 11 FPM 得 到 的 不 同 冲击 动量 下 点 C.D.E 的 位 移 包 络 图 
图 11 以 第 一 指标 判断 动力 失 稳 ,由 图 可 知 : Fig.11 Displacement envelope of point C,D and Ek 


under different impact momentum simulated by FPM 


1) 当 冲 击 动量 达到 120 kg * m/s 时 3 个 节点 位 移 包 
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图 12 以 第 2 指标 判断 动力 失 稳 。 图 中 给 出 不 
同 动量 下 点 CD、E、P 的 最 大 竖 向 位 移 对 比 。 对 于 
点 CC 点 DD 和 点 , 当 冲 击 动量 为 120 kg . m/s 时 节 
点 竖 向 位 移 最 大 值 大 幅 增 加 ,节点 发 生 失 稳 。 但 对 
于 点 下 ,试验 所 得 位 移 在 动量 为 20 kg . m/s 就 已 失 


4.36 5 62 8.72 10 20 120 
动量 /(kg*m*s') 
(a) FPM 
0 
人 

2D ©. 一 
E 
: 
长 -20 
哈 
浪 
要 
=30 eC 
i 一 e 一 点 D 
Pp 一 人 一 点 E 


4.36 5 62 372 10 20 120 
动量 /kg .ms ) 
(b) 试验 
12 不 同 冲击 动量 下 点 C.D、E、F 最 大 坚 向 位 移 对 比 
Fig.12 Comparison of maximum vertical displacement of 


point C,D,E and F under different impact momentum 
4.4 能 量 传递 过 程 分 析 


单 层 柱 面 网 党 的 能 量 传递 分 为 2 个 研究 区 段 ， 
如 图 13 所 示 , 以 工 况 7 为 例 说 明 能 量 传递 过 程 。 如 
图 13(a) 所 示 , 工 况 7 冲击 结束 时 间 为 93. 06 ms ,发 
生 冲 击 后 ,冲击 区 杆 件 2 和 杆 件 3 率先 产生 动力 响 
应 。 随 后 依据 距 撞 击 点 距离 , 其余 杆 件 依 次 发 生 响 
应 ,但 整体 响应 时 间 不 超过 5 ms, 远 小 于 冲击 结束 时 
间 , 其 中 以 非 冲击 区 的 杆 件 1 动力 啊 应 最 为 明显 ,说 
明 在 冲击 结束 前 ,网 壳 结 构 整 体 已 经 产生 了 明显 的 
动力 响应 。 同 时 在 这 一 阶段 ,能 量 更 易 朝 着 靠近 结 
构 边缘 的 杆 件 1 和 杆 件 2 传递 ,而 同 为 冲击 区 的 杆 
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件 3 的 应 变 能 较 杆 件 2 要 小 。 这 是 因为 撞击 点 左 侧 
结构 的 纵向 长 度 小 于 其 右 侧 结构 , 故 其 刚度 更 大 , 受 
力 时 所 分 配 的 力 更 多 ,应 变 能 更 大 。 
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Fig.13 Changes of strain energy of members 
图 13(b) 则 说 明 能 量 更 易 沿 着 网 壳 斜 杆 传 递 。 
第 2 区 段 杆 件 的 应 变 能 幅 值 在 0.361 ~4.86J 之 间 ， 
显著 大 于 第 1 区 段 杆 件 应 变 能 幅 值 范围 , 杆 件 6、8 
和 9 产生 的 能 量 均 大 于 杆 件 1。 距 离 支 座 较 近 的 杆 
件 6 和 处 于 冲击 区 的 杆 件 8 在 随后 的 振动 阶段 的 应 
变 能 出 现 多 个 峰值 ,说 明 单 层 柱 面 网 诗 受 冲击 丛 载 
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作用 时 ,冲击 区 和 支 座 为 受 力 的 薄弱 部 位 ,其 稳定 性 
较 差 。 图 13(c) 可 印证 上 述 结论 ,在 跨度 方向 、 支 座 
周边 和 冲击 区 周边 的 杆 件 能 量 较 大 ,同时 斜 杆 所 产 
生 的 应 变 能 要 大 于 纵向 杆 件 ,这 与 基于 力 流 理论 的 
单 层 球面 网 壳 结构 的 传 力 特点 具有 相似 性 ” 。 


4.5 误差 分 析 


上 述 对 比 中 ,有 限 质 点 法 能 够 反映 真实 情况 下 
网 壳 结 构 受 冲击 时 的 动态 响应 ,并 能 从 中 推导 出 与 
实际 情况 相符 的 结论 ,但 在 部 分 数值 ,尤其 是 位 移 部 
分 ,有 限 质 点 法 和 试验 的 差 值 较 大 ,原因 分 析 如 下 。 

1) 有 限 质 点 法 未 能 考虑 网 充 结 构 焊 接 节 点 的 残 
余 应 力 ,尚未 建立 网 壳 结 构 节 点 损伤 模型 ,致使 有 限 
质点 法 网 壳 模 型 有 更 好 抗 冲击 性 ,所 得 的 节点 竖 向 
位 移 和 杆 件 应 变 要 小 于 试验 结果 。 

2) 有 限 质 点 法 采用 的 对 称 罚 函 数 碰撞 模型 简化 
了 冲击 面 几 何 形状 、 粗 糙 度 .接触 面积 的 影响 ,以 冲 
击 物 和 被 冲击 物 发 生 穿 透 计 算 冲 击 力 。 

3) 有 限 质 点 法 采用 点 值 描述 ,将 杆 件 的 质量 全 
部 集中 在 质点 上 ,单元 不 存在 弯曲 变形 ,弱化 了 网 壳 
结构 节点 受 杆 件 的 牵连 作用 。 


5 结 论 


1) 本 研究 根据 能 量 守恒 定律 提出 了 基于 有 限 质 
点 法 的 结构 碰撞 刚度 系数 确定 方法 。 利 用 以 对 称 罚 
函数 为 基础 的 双 球 碰撞 模型 建立 有 限 质 点 法 冲击 模 
型 , 编 入 适用 于 解决 金属 冲击 问题 的 Cowper-Sym- 
onds 本 构 模 型 。 

2) 对 比 不 同 工 况 下 单 层 柱 面 网 壳 定 点 受 冲击 作 
用 的 动态 响应 ,结果 表明 :中 从 冲击 力 幅 值 来 看 ,有 
限 质点 法 获取 的 冲击 力 幅 值 与 试验 吻合 度 较 高 , 冲 
击 力 幅 值 的 误差 范围 为 -20.29% ~ 14.98% ;从 冲 
击 答 载 形 式 来 看 ,有 限 质 点 法 时 程 曲线 近似 半 正 疏 
曲线 ,试验 所 得 时 程 曲 线 近 似 三 角形 曲线 ,经 验算 冲 
击 荷 载 形式 不 会 成 为 影响 结构 动态 响应 的 主要 因 
素 。 忆 从 动态 失 稳 分 析 来 看 ,结构 的 整体 失 稳 动量 
为 120 kg ' m/s, 靠 近 支 座 的 节点 玉 受 文 座 约束 的 影 
响 ,其 竖 向 位 移 要 大 于 其 他 3 个 节点 。 同 时 , 工 况 4 
的 动力 响应 要 小 于 冲击 动量 更 小 的 工 况 1, 表 明 冲 
击 造 成 的 结构 动态 响应 在 一 定 程度 上 由 冲击 速度 控 
制 。@ 从 能 量 传递 过 程 来 看 ,冲击 产生 的 应 变 能 更 
易 沿 着 网 壳 斜 杆 传递 ,并且 冲 击 区 杆 件 和 刚度 较 大 
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的 杆 件 如 路 度 方 向 ` 支 座 处 也 成 为 主要 的 受 力 构 件 。 

3) 有 限 质点 法 能 模拟 出 结构 受 冲击 时 的 动态 响 
应 ,所 得 冲击 力 、 节 点 位 移 、 应 变 等 参数 的 时 程 变化 
规律 与 试验 一 致 ,可 以 从 定性 和 定量 2 个 方面 判断 
结构 受 冲击 的 薄弱 部 分 ,为 提高 相关 结构 抗 冲击 性 
能 提供 理论 依据 。 
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